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The probability of error in the detection of RF pulses embedded in a com- 
bination of Gaussian noise, time crosstalk from the tails of two neighboring 
pulses, and frequency crosstalk from an adjacent channel, is calculated. 

It is shown that for a given probability of error it is possible to maximize 
the pulse repetition frequency and simultaneously to minimize the channel 
spacing and signal-to-thermal noise by operating the system at a signal-to- 
thermal noise level close to the level of the combined time and frequency 
crosstalk. 

I. INTRODUCTION 

Consider many PCM messages occupying adjacent frequency bands 
in the same transmission medium, as, for example, in the proposed long 
distance waveguide communication system. 1 At some point the messages 
must be separated and read; these operations are performed by the re- 
ceivers. Each receiver will be considered to consist of a filter and an en- 
velope detector that takes periodic instantaneous samples and decides if 
the level of the signal is above or below a threshold. 

Suppose that at a certain sampling time there is no pulse to be de- 
tected. Nevertheless, the received signal will be composed of the summa- 
tion of three types of interference: time crosstalk or intersymbol inter- 
ference, frequency crosstalk, and noise.* 

Time crosstalk is measured by the envelope of the message at the 
sampling time in the absence of other messages and noise; it is due to 
the trailing and leading edges of the other pulses that make the message. 
It is known that if the sampling is instantaneous the time crosstalk can 
be reduced to zero by proper choice of filters and input signal, 2 but in a 



Throughout this paper we understand "noise" to be thermal noise. 
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real system the sampling time is not zero, and consequently the inter- 
symbol interference varies during that time. The actual description of 
how this varying crosstalk influences the detected signal is a very com- 
plicated problem that involves a detailed knowledge, not only of the 
input pulses and transfer characteristics of transmitters and receivers, 
but also of the detector. We by-pass this problem by assuming conserva- 
tively a fictitious system that indeed has instantaneous sampling, but 
with the time crosstalk being the maximum value achieved by the time 
crosstalk in the real system during the finite sampling time. 

Frequency crosstalk is measured by the envelope at the sampling 
time in the absence of the wanted message and the noise; it is due to 
the fact that the other messages have spectrums that overlap with the 
transfer characteristic of the receiving filter of the channel under con- 
sideration. 

Finally, noise is measured by the envelope at the sampling time in 
the absence of all the messages; it comes essentially from the first ampli- 
fier in the receiver. 

If the envelope of the three interferences is bigger than the slicing 
level, the detector decides that a pulse exists in that time slot, and an 
error is made. Similarly, the detector makes another error if a pulse 
should be detected but is shadowed by the interferences in such a way 
that the envelope of the received signal is smaller than the slicing level. 

It is the purpose of this paper first to determine the relationship be- 
tween the amplitudes of the wanted signal, time crosstalk, frequency 
crosstalk, noise, and slicing level; and second to establish in some sense 
the most efficient design of a system for a given probability of error. 

II. DENSITY DISTRIBUTION OF SIGNAL, TIME CROSSTALK, FREQUENCY 
CROSSTALK, AND GAUSSIAN NOISE 

Simplifying assumptions: 

(a) Time crosstalk is represented by the sum of two sine waves of 
the same amplitude and arbitrary phases. The implications are: First, 
only the trailing edge of the preceding pulse and the leading edge of the 
following one are important. Second, each received pulse is symmetrical. 
This is rigorously true if the input pulse is symmetrical and the system 
has no phase distortion. Third, the phases of the pulses are uncorrelated, 
which is true if the pulses have passed through several partially regen- 
erative repeaters. 

(b) Frequency crosstalk is represented by a sine wave of arbitrary 
phase. The implications are: First, only one neighboring channel feeds 
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non-negligible power into the wanted channel. This is shown to be a 
reasonable assumption in Appendix A. Second, the pulses in different 
channels are synchronized. If they were not, the amplitude of the fre- 
quency crosstalk would vary between the two extreme values that can 
be obtained with the best and the worst interleaving of pulses. 

(c) The noise is assumed to be Gaussian. 

(d) The detector measures instantaneously if the envelope of. the 
received signal is above or below a threshold. This is probably the crudest 
approximation, because in a real system the detector is not ideal and, 
what is even worse, the signal passes through repeaters with only partial 
regeneration. 

The vector representing the signal to be detected is 

*S' = A + p r c' ' + p T c' " + pfc' 3 + Gaussian noise, (1) 

where A is the amplitude of the RF of the wanted pulse; its value is 
one if there is a pulse to be detected, and zero if there is no pulse; its 
phase is taken as reference. The second and third term represent the 
time crosstalk; they arc vectors of the same modulus p T , but arbitrary 
phases 81 and Q* . The fourth term represents the frequency crosstalk of 
modulus pf and arbitrary phase 6 S . Each one of these three last vectors, 
being originated from binary pulses, has a 50-50 chance of being present 
or not. The bivariate density distribution,* calculated in Appendix B, 
(52) is 

+ c~»' W h (^') + 2,- ( >^ 2/2 '7 (^) h (^') (2) 

where x and y are the coordinates of the terminus of S, the vector repre- 
senting the signal to be detected; /• = \/(x — A ) 2 + j/ 2 ; a 1 is the mean 
noise power; and I is the modified Bessel function of first kind of order 
zero. The density distribution (2) is only valid for the tail of the dis- 
tribution ; that is, 

Pt ( 

«r. (3) 

Pr 



For tutorial background see, for example, Bennett. 3 



924 



THE BELL SYSTEM TECHNICAL JOURNAL, MAY 1961 



It is possible to interpret the meaning of each term in (2). The first 
one is the contribution to p(x,y) when only noise is present; on the aver- 
age, this combination happens once each eight detections. The second 
term is the contribution when noise and only one of the two time cross- 
talk tails are present; on the average, this combination occurs once each 
four detections. The third term is the contribution when noise and both 
time crosstalk tails are present; on the average, this combination occurs 
once each eight detections. The fourth term is the contribution when 
noise and frequency crosstalk are present ; on the average, this combina- 
tion occurs once each eight detections. The fifth term is the contribution 
when noise, frequency crosstalk, and one time crosstalk tail are present; 
on the average, this combination occurs once each four detections. The 
sixth term is the contribution when noise, the two time crosstalk, and 
frequency crosstalk are present; on the average, this combination hap- 
pens once each eight detections. 

If there is a pulse to be detected (pulse on), A is equal to one and the 
density distribution (2) is 

p(x,y) = pi(x,y) (A = 1). (4) 

If there is no pulse to be detected (pulse off), A is zero and the density 
distribution is 



(A = 0). 



(5) 



p(x,y) = pi(x,y) 

Both functions, Pi(x,y) and Pi{x,ij), schematically plotted as Figs. 1(a) 
and 1(b), have the same bell shape and circular symmetry around their 
respective axes located at x = 1, y = 0, and at x = y = 0. 




p 2 (*.y) 



(a) 




Fig. 1 — Density distribution for (a) pulse on, (b) pulse off. 
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III. PROBABILITY OF ERROR 

In general, the volume defined by 

P = [ p(x,y) dx dy (6) 

measures the probability that the signal S be a vector originating at 
the origin of coordinates and terminating at any point within the area 
of integration a. 

The quantity (1 — P) measures the probability that £ is a vector 
with the terminus outside the area a. The detector decides whether the 
terminus is inside or outside of a. 

Suppose that the signal free of interference A has its terminus out- 
side of the area a; then if the received signal S is also outside of a the 
detector makes a correct decision, but if S falls inside of a the detector 
makes an error. Since the probability of finding S inside of a is given 
by P, this integral measures the probability of error and (1 — P) meas- 
ures the probability of making a correct decision. 

The detector we use is one capable of deciding if the envelope of the 
received signal is bigger or smaller than a threshold p . 

The probability of error in the "on pulse" condition, Fig. 1(a), is the 
probability that \ S \ < p : 

Pi = [ Pi(x,y) dx dy , (7) 

where pi(x,y) is derived from (2) by setting A = 1, and ai is the circle 
of radius p and center at the origin of coordinates. The integration per- 
formed in Appendix C yields (70): 



K, 



A = 



d -po)H 

2<j2 J[j +26- p r' /2 '7o(, ' 



L6ir 



n I Pt 



(8) 



+ (' 7 h I Pt j — ) + e * JoIpf g 

where 7 is the modified Bessel function of the first kind of order zero 
and K is the modified Bessel function of the second kind of order zero. 
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The probability of error in the "off pulse" condition, Fig. 1(b), is the 
probability that \ S\ > po : 

P 2 = f p 2 (x,y)dxdy, (9) 

where f*{x,y) is obtained from (2) by making ^4=0, and a 2 is the 
surface outside a\ . The integration performed in Appendix C yields 

(77): 

-p 2/2 ff 2 



P,= 6 - 



[l + 2<fV*»J. (e^j + e+r'"*!* (^) 



+ e-" w / ( *■=? ) + 26- l '^' J W / B ( ^- u ) /, ( ^ u ) (10) 



+ -Up T +,y)*l1**ji /prPo\ ^ /PFP0\ 



with Jo and Ko being the modified Bcssel functions of the first and 
second kinds. 

The six terms appearing in expressions (8) and (10) have the same 
physical interpretation as that given for the six terms appearing in (2). 

Since the "on" and "off" pulses are equally likely, the probability of 
error of the message is 

p = kp, + P 2 ). (11) 

The probability of error of the message P can be calculated for any 
combination of time and frequency crosstalk p T and p F , but it is possible 
to relate these two values by demanding that, according to some rule, 
both are equally damaging to the system. The rule we adopt is given 
by the following equations: 

Pi(pf = 0, po = 0.5) = Pi(pt = 0, po = 0.5); 
P 2 (pf = 0, po = 0.5) = P 2 (pr = 0, po = 0.5). 

For any signal-to-noise level and a slicing level equal to half the pulse 
amplitude (po = 0.5), the probability of error in the "on" or "off" 
pulse condition due to noise and only time crosstalk is equal to that 
due to noise and only frequency crosstalk. 

Substituting (8) and (10) in equation (12), we get 

2e->r 2/2 '7„ (V) + e- 2 ^ 2/ 'V fe) - 1 = 2e- p » w h fe) . (13) 
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Frequency crosstalk p F has been plotted against time crosstalk p T , 
for different signal-to-noise levels l/\/2o- in Fig. 2. 
A line defined by the following equation 

20 log - = 20 log — + 3 

Pt Pf 

has been included in the same figure (dotted line) for comparison pur- 
poses. Either from (13; or from Fig. 2 it can be deduced that for 



Pt 

0^2 



«1, 



Pf 



~V2 



Pt 



(14) 



and for 



Pr 



»1, 



Pf = -Pt- 



(15) 
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Fig. 2 — Equally damaging time and frequency crosstalk. 
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Fig. 3 — Probability of error of a message in the presence of crosstalk and noise. 
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For a given signal-to-noise ratio, if the normalized time crosstalk in- 
tensity per tail, pr , is small compared to the normalized mean noise 
power 2a 2 , two time-crosstalking tails introduce by themselves as many 
errors as does one frequency crosstalk 3 db above the level of each tail. 
But if p T » 2o 2 the time-crosstalking tails introduce as many errors as 
does one frequency crosstalk G db above the level of each tail. 

Tor each set of values a, p T , and p F that satisfies (13) we calculate 
from (11) the optimum slicing level po that minimizes the probability 
of error, and P m i n , the value of that minimum. Fig. 3 contains this 
information. The probability of error, P m in , is plotted as a function of 
signal-to-noise level for different values, p T , of time crosstalk per tail. 
Each set of pairs of numbers on these curves indicates the local optimum 
slicing level p and the frequency crosstalk p F . 

The dashed line (no crosstalk) almost coincides with that derived by 
Bennett. 3 The small difference stems from the fact that Bennett calcu- 
lates the probability of error of the message for equal contributions of 
errors from the "on" and "off" pulse condition, while we calculate the 
minimum probability of error of the message. 

IV. OPTIMUM DESIGN REGION 

Suppose that we want to design a system with a given probability of 
error. Is there only one combination of values of crosstalk and signal-to- 
noise capable of satisfying the demanded probability of error? The 
answer is no. In Fig. 7 the given probability of error will be an ordinate 
obtainable with an infinite number of combinations of signal-to-noise 
level and crosstalk. We will develop two criteria for making a reasonable 
choice, and for that purpose we need some intermediate steps. 

As a first step we redraw the part of Fig. 3 for low probability of error 
in Fig. 4, using time crosstalk per tail as the abscissa, signal-to-noise 
as the ordinate, and probability of error as parameter. The frequency 
crosstalk and optimum slicing level change slightly from point to point, 
but their exact values have not been written down. 

As a second step Ave derive Fig. 5 from Figs. 6 and 7, which, together 
with Fig. 8, a sheet of definition of symbols, have been taken from the 
companion paper. We shall see later how the derivation takes place, 
but. first let us get acquainted with Figs. 6 and 7. In both these figures, 
the ordinates are proportional to time spacing between successive 
pulses, t, times frequency spacing between adjacent channels, \f\ — f 2 \ . 
The smaller this product the better, because time and frequency oc- 
cupancies arc proportional to the product. The different coefficients of 
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Fig. 4 — Reproduction of a section of Fig. 3, using new coordinates. 



proportionality in the figures have to do with the input pulse width 2T 
and sampling time 2T r . The abscissas measure the ratio between band- 
width of the sending filter, 2Fi , and the bandwidth of the receiving fil- 
ter, 2/^2 • Eack figure contains three sets of curves, each corresponding 
to different time crosstalk per tail and different frequency crosstalk. 
Finally, the curves in each set correspond to different combinations of 
transfer characteristics of the transmitting and receiving filters. 

The upper and lower dashed lines in Fig. 5 are applicable to systems 
with sending and receiving filters, each approximately maximally flat 
(three cavities); they have been derived from the dotted lines in Figs. 
6 and 7, respectively. The upper and lower solid lines in Fig. 5 are ap- 
plicable to systems with Gaussian sending filter and receiving filter 
approximately maximally flat (three cavities); they have been derived 
from the full lines in Figs, b* and 7, respectively. The ordinates in Fig. 5 
are the ordinates of the minimums of Figs. 6 and 7, and the abscissas 
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iii Fig. 5 are the different time crosstalks per tail corresponding to each 
set of curves in Figs. (> and 7. 

It is important, to bear in mind that the ordinates of Fig. 5 are pro- 
portional to the minimum lime spacing, t, times channel frequency 
spacing, |/i — /l>| , which corresponds to maximum rate of information 
transmission. 

As a third step we compare Fig. 4 with Fig. 5. For the same value of 
the abscissa both figures have ordinates that measure properties of the 
system we want to be as small as possible, but, since the slopes in the 
two figures are of different sign, a system operating at high time cross- 
talk per tail (small abscissa) will have (Fig. 5) a desirable low value 
t|/i — f 2 \ but a large and unwanted signal-to-noise level. Conversely, a 
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Fig. — r | /, - /. | curves lovr/2T = 1.5 + T r /2T. 

system operating at low time crosstalk per tail will have an undesirably 
large r|/i — /2I and a wanted low signal-to-noise level. This suggests 
the existence of an intermediate optimum, and the question now is what 
function we want to minimize. 

The answer is elusive, because what we really want is to minimize 
the price of a system that handles a certain rate of information with a 
given probability of error. That cost must be a function of time spacing, 
channel spacing, signal-to-noise ratio, and perhaps other variables. We 
don't know that function — at least not now — and because of lack of 
better knowledge we propose the minimization of two simple functions 
in which the signal-to-noisc ratio is weighted differently: 
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where 20 log (l/V^tr) is the signal-to-noise in db and \/y/2o is the 
ratio of rms signal and rms noise. 

The different weighting functions were selected in order to introduce 
some, idea about the influence of distance between successive repeaters. 
Since the amplitude of the received signal decays exponentially with the 
distance, between terminals, for fixed transmitter and receiver G'i de- 
creases linearly with distance and G 2 decreases . xponentially with 
distance. 

Functions G\ and G 2 , obtained by multiplying the ordinates of each 
curve in Fig. 5 by the properly weighted ordinates of Fig. 4, have been 
plotted in Figs. 9, 10, 11, and 12. Each figure contains two sets of curves, 
and in each set the three curves exhibit minimums individualized by 
the coordinates probability of error and time crosstalk per tail. Those 
coordinates identify three points of an optimization curve that could 
be plotted in Fig. 3. For clarity, part of Fig. 3 has been reproduced in Fig. 
13, omitting the detailed information on frequency crosstalk and opti- 
mum slicing level. In Fig. 13 we have plotted the lines joining each set of 
three points rather than the points themselves. Since there are eight 
sets of curves in Figs. 9 through 12, we get eight lines of optimum design 
in Fig. 13. Four of them correspond to the minimization of G\ (signal- 
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Fig. 9 — Minimization of G\ , transmitting and receiving filters approximately 
maximally flat (three cavities). 

to-noise in db) and arc clustered close to the line defined by the param- 
eter time crosstalk per tail 24 db; the other four lines of optimum design 
correspond to the minimization of G 2 (rms signal to rms noise.) and are 
close to the line defined by the parameter time crosstalk per tail 20 db. 
In each cluster, the two solid lines are optimization curves for two 
systems, both with maximally flat (three cavities) transmitting and 
receiving filters but with different input pulse width 2T and sampling 
time 27' r ; the two dashed lines are optimization curves for two systems 
both with Gaussian transmitting filters and maximally flat (three 
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Fig. 10 — Minimization of G t , transmitting filter Gaussian, receiving filter 
approximately maximally flat (three cavities). 

cavities) receiving filters, but with different input pulse width 2T and 
sampling time 2T r . 

In spite of the different dependence of signal-to- noise in G\ and G 2 , 
the different shapes of transmitting and receiving transfer characteristics, 
and the different input pulse widths and sampling times, all curves of 
optimum design are rather close to each other, and they are essentially 
located in the region where rms noise and rms crosstalk are comparable. 

The optimum design lines are in general slightly steeper than the 
constant time crosstalk per tail lines, and, in particular, the two ex- 
tremes of each of these eight design lines correspond to a change of 100 
in the probability of error, around 1.5 db in signal-to-noise and only a 
few tenths of a db in time crosstalk per tail. This means that, once an 
optimum system has been built, the probability of error can be changed 
substantially by modifying only the signal-to-noise ratio (which is easy 
to do) and, in spite of this change, the system will remain close to the 
optimum design. 
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Fig;. 11 - Minimization of G-< , transmitting filter Gaussian, receiving filter 
approximately maximally flat (three cavities). 



V. CONCLUSION'S 

The probability of error in the envelope detection of an RF signal 
embedded in Gaussian noise, time crosstalk from two neighboring pulses, 
and frequency crosstalk from an adjacent channel has been calculated 
and plotted in Fig. 3. 

Also, two kinds of optimum operating conditions have been postulated 
which yield the results shown in Fig. 13. These conditions allow one to 
design a system in such a way that some minimization of time spacing 
between successive pulses, frequency spacing between adjacent channels, 
and signal-lo-noise ratio is achieved. 

An example of design is this: Suppose wc want an optimally designed 
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Fig. 12 — Minimization of G 2 , transmitting and receiving filters approximately 
maximally flat (three cavities). 



system that has a probability of error of 10~ . We don't know which of 
the two criteria of optimization developed in this paper is closer to 
reality, and, because of that lack of knowledge, we adopt the middle 
of the road for the example. In Fig. 13, the ordinate 10~ and the middle 
of the optimum design region establish that the system should have a 
signal-to-noise level of about 20.6 db and a time crosstalk per tail of 25 
db. This last datum is enough to enter in the companion paper and to 
complete the design of the system. 
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Fig. 13 — Optimum design region: solid curves — transmitting and receiving 
filters approximately maximally flat; dashed curves — transmitting filter Gaus- 
sian, receiving filter approximately maximally flat. 



APPENDIX A 



Crosstalk Between Adjacent Frequency Channels 

We want to determine the frequency crosstalk between adjacent 
channels in order to find what arrangement of filters is the most favorable. 
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Fig. 14 — Arrangements of transmitting and receiving niters in a system. 

The systems we shall deal with, shown in Figs. 14(a) and 14(b), differ 
in the order in which the bands arc dropped. Each system consists of 
many transmitting and receiving filters, of which only three transmitting 
filters — To , Ti , and t 2 — and three receiving filters — Ro , Ri , and 
r 2 — arc drawn, because we assume that the crosstalk in a receiver (rj), 
comes essentially from the immediately neighboring channels. 



Fig. 15 — Transmitting or receiving filter. 

For simplicity, we assume that, except for the frequency at which 
they are tuned, all the filters are similar to that shown in Fig. 15. They 
are constant-resistance, symmetrical, and reciprocal, and the transfer 
functions between terminals are given by the scattering matrix 

»J11 »^12 »->l3 Ol4 



s = 



021 022 '->23 
031 032 033 
0,11 042 ^'43 


024 




iy/\ 


Y 






Y 







?Vl - Y' 2 


Y 




i'Vl - Y 2 




Y 


iy/\ - 



Y 2 





- Y 2 





(16) 
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Now we can calculate the maximum intensity of the crosstalk C mq 
between m and q in Fig. 14(a) due to a pulse entering at m. 

Assuming for simplicity that the system is phase equalized and that 
the input pulse has a (sin x)/x shape (rectangular spectrum), the 
maximum intensity of the crosstalk is given by 



C, 



= K [ G»Y n y/\ - Y? Vl - Y-? Vl - YjYidf, (17) 



where K is a constant of proportionality, G is the rectangular spectrum 
of the input pulse centered at/ , and Y \/l — Y{ 2 \/l — Y Y\ derived 
from (16), and Fig. 14(a) is the transfer function between m and q. 
The subindices and 1 refer to the center frequencies / and /i of each 
scattering coefficient. 

Following arguments similar to the preceding one, the maximum 
crosstalk intensities between n and Q in Fig. 14(a) and between s and v 
and between u and v in Fig. 14(b), are 

C NQ = K I G»Y, Vl - Y<? F x df, (18) 

Jo 

C sv = K I G Q Y Vl - 5'V 2 (1 - Y 2 i )Y 1 df, (19) 

Jo 

C uv = K [ G 2 Y-> y/\ - Y? Y 1 df. (20) 

Jo 

The factors involved in each integrand of (17) through (20) have 
been plotted in Figs. 16(a), (b), (c), and (d). The integrands of (17) 
through (20) are obtained by multiplying the curves in Fig. 16(a) 
through Fig. 16(d) respectively. The results which happen to be the 
output spectra are plotted in Fig. 17(a) through Fig. 17(d). 

The integration of these curves with respect to frequency, that is, the 
areas between the curves and the frequency axes, are, because of (17) 
through (20), proportional to the maximum intensities of the crosstalks. 

Comparing these areas we deduce 

(- MQ = C UV , (21) 

(-■KQ = CsV , \^-) 

C MQ «C* Q , (23) 

C uv «C av . (24) 

The first, two equations show that the total crosstalk in the system of 
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Fig. 10 — Factors in volved i n integrands of (17) through (20); G H = input 
signal spectra; }',„ , vl — IV = scattering coefficients. 

Fig. 14(a) is the .same as the total crosstalk in the system of Fig. 14(b). 
Furthermore, from (23) and (24) we deduce that the total crosstalk in 
either system comes from the superposition of two signals, of which one 
is negligible compared to the other. 
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Fig. 17 — Spectra of crosstalk signals. 



APPENDIX B 



Bivariate Density Distribution 

An on-off pulse embedded in unwanted crosstalk and noise is repre- 
sented vectorially by 



S = A -\- pre 1 ' 1 + Pre' 9 ' + Pfc' 83 + Gaussian noise. 



(25) 



The amplitude A of the RF pulse is unity if the pulse is on, and zero 
if the pulse is off; the phase of this vector is selected zero, as reference. 
Time crosstalk is represented by two vectors of the same modulus p T 
and arbitrary phases 0, and 2 . Frequency crosstalk is represented by a 
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vector of modulus p F and arbitrary phase 3 . Each one of these three 
vectors, since it originated from binary pulses, has equal probability of 
being present or not. The phases d\ , 6 2 , and 8 3 have a constant probability 
of acquiring any value between zero and 2r. 

We want, to calculate the density distribution of <S, and we know the 
density distribution of each one of its five uncorrelated terms: 

Pl (x,y) = S(x - A)S(y), (26) 

»(**) = »<**) - * {V ^±^ + 5(V - 7 T ? ~ Pr) , (27) 
■±7r V •' + if' 47r Pt 

, ... u , = KV^±£l + '(^TZ^f) , ( , 8) 

47rV-r- + y- 47r Pf 

p*(x,y) = 1 e -c^ W> (29) 

27TCT" 

where 6(2) is the Dirac delta function and a is the variance, which in 
this particular problem measures the mean noise power. 

It is known that the distribution p(x,y) of the sum »S' of independent 
terms is equal to the inverse transform of the product of the double 
Fourier transform of the density distribution of each term of the sum. 3 

The double Fourier transform of a distribution p„(x,y) is, by defini- 
tion 



'■ - // 



e n " x+,rj) p n (x,y)dxdy. (30) 



Replacing p„(x,y) by (20) through (29) and integrating, 

C x = e*\ (31) 

(\ = C, = J[l + MptVFT^)] , (32) 

r, = f[i + j„UVF+^-)]. (33) 



-(«* 



Ct - 6-™*™ t (34) 

where Jo is the Bcsscl function of first order and kind. 

The inverse transform of the product of these functions is the distribu- 
tion p(x,y) of the signal S we were looking for: 

v(x,v) = JL (/>"«*-' >+».-<.■»>«*+.»[, + UerVF+7)]' ,. , 

•">-7T" J J (35) 

— OC 

■[1 +./o(p,VFT7 2 )]^^- 
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The distribution p(x,y) is the summation of several double integrals 
to be evaluated. The most general of them is 

_t (36) 

■MpzV?- + r)Jo{pW? + r,-) d£dr,. 
Replacing each Bessel function by an integral expression, 5 

lit 00 

W = — L fff dadddy [[ c-'^-^-Ol^a-p^osff-p^osyl 

(Iir) 6 JJJ JJ (37) 

d^dr,. 





— »ij(tf— Pi sin a—p-i sin 0— p 3 sin 7) — (o 2 /2) ($ 2 + ij 2 ) 



The two integrations from — oo to °o arc known Fourier transforms, and 
(37) becomes 

W= l 



(27TO-) 2 

trr • (38) 

• / / / «' (j:_A_P1 cos a_p - cos P~ Pi cos 7> (""Pl 8,n a ~''2 am 0-p 3 sin 7) 2 ) /2ff 2 j j« , 


By changing variables, 

.r — -1 = r cos ys, 

(39) 
y — r sm y 

the exponent can be rearranged: 

-(r 2 +Pl 2 +P2 2 +p 3 2 )/2,r2 

W = - .^-^ 



III Irpi cos (a— ?)+rp« cost (fl— <p)+rp 3 con (y— if) ( Af)) 

JJJ -P1P2 cos (a-0)-pip 3 cos («— 7)— P2P3 cos (0-7>]/<r 2 j ,„ . 

We start integrating with respect to a. The integral to be solved is 
essentially 

TT/- __ / p Pli r ros (a— V)— P2 cs (a— 0)— p 3 cos (a— 7)] /a 2 j / , * \ 

JO 
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The exact result is 



W a = 2tIo[ h -\ 



pi 



[r — pi cos (cp — j8) - p 3 cos (<p — y)] 2 

4- [p 2 sin (<p — /3) + ps sin (<? - 7)]' 



(42) 



but if we carry this expression to (40) the integration with respect to 
/3 and 7 becomes extremely complicated. 

A substantial simplification can be obtained if we consider first that 
we are interested only in the tails of the distributions, and consequently 



r »< 



pi 
p-i 

\p-i 



(43) 



Second, for 



rpi 
a 2 

rpz 

^2 



»1, 



(44) 



which is the only nontrivial case, the main contribution to the triple 
integral (40) comes from values of the integrating variables 



r = <P- 



(4; 



Because of (43) and (45), expression (41) can be reduced to 
W a ^ e -' 1 "' ,,+ " ) ' g ' 1 C e (Pir/ff2> c ° 8 {a ~ v) dip 
and, after performing the integration, 



(46) 



(47) 



The reader may also derive this result from (42), (43), and (45). 
Substituting this result of the integration on a, in (40), 



ERRORS IN DETECTION' OF RF PULSES 947 

/ ,(r2+p 1 2+p 2 2+p 3 2+2p 1 pj+2p l p 3 )/2 1 r2 / v 

"' = e g-i /.r) 

27RT- \ (T- / 

(48) 

[rp 2 cos (0-^+rp;, cos (T-p)-p-2Pa <--°s (P-?)]/* 2 Jg J 

II 

Now we perform the integration on /3. The integral to be solved is 
essentially 

W, = f^ e Wff2)tr COB ( ^- p3 coa (P ~ y)] d3. (49) 

Following the same reasoning used to integrate W a , in (41), the 
approximate result is 

W 9 » J^-"-'" 7o (^) . (50) 

After substituting in (48) and performing the integration on y, W is 

W = £ e -^+>2+,3> W /() ^ /o ^ /() M (51) 

Substituting this generic result in (35), the density distribution of the 
signal S is obtained : 

+ e -" 2/2 ' 2 /„ (^) + 26- 1 '**"™* h (^*) /o (^) (52) 

APPENDIX C 

Evaluation of Probabilities of Error 
Case 1: "Pulse On" 

We want to evaluate the integral 

P x = / px (.<:,//) dx dij. (53) 

J "i 

where a\ is a circle of radius p (l ^ ^ with center at the origin of co- 
ordinates, and pi(x,y) is obtained from (52) by making 

A = 1 (54) 
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ill the expression 

r = V(x - AY + y*-. (55) 

Adopting the following change of variables: 

X = p cos \I/, 

(56) 
y = p sin \p, 

Pi becomes 

I\ = I I pi(p cos \p, p sin \f/) dp dip. (57) 

The most general term of the integration is proportional to 

^-if'-^feWSW?)'** (58) 

where 

r = Vl + p 2 - 2p cos 4f. (59) 

We simplif}' the integrand. Notice first that, since 

Po £* I, (60) 

we deduce, from (59), 

independently of yp; second, the range of interest for a is 

<r«po. (61) 

Therefore, the exponent in (58) is 

f-, » I- (62) 

Bccause of this inequality and because, for a small variation of r, the 
exponential in (58) varies much faster than the modified Bessel func- 
tions, most of the contribution to the integral comes from values of the 
variable close to those that minimize r, 

* = o, 

P = Po , 
and (58) becomes 

Ux £* / ( Pl L^A°) / ( P2 Lr ») /o ( p:t - 1 " ») D> (63) 

D _ J* /"° e- (1+p2 - 2p cos < )/2 * 2 P <ty dp. (64) 



where 
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Integrating with respect to rf/, 



1) 



= 2tt J e in I -., I p dp. 



(65) 



Since the exponential varies much faster than the rest of the integrand, 
most of the contribution to the integral comes from p = po , and, be- 
cause of (60) and (61), hip/a 2 ) can be replaced by its asymptotic 
expansion . Consequently, 



rPti 

D ^ \Z2irpo o I e 



-(l-p)2/2 ff 2 



<h. 



Integrating, 



D = y/2^o.-^— C - (1 - p0,2/2ff2 - 



1 — pn 
for compactness, this can be rewritten 

^ 2 K r d - po) 2 i 



D 



(66) 



(67) 



:68: 



where 7f is the modified Bessel function of the second kind. 

Substituting (68) in (63), the general term of the integration (57) is 



Pn 



and the probability of error for the "pulse on" condition is 

I 2<r 2 r. . ., - P -.2/2.72 r / 1 — p<>\ 

1 — =* 1 + 2< 7 in I Pr — ^ — I 

1 — Pn 

■ e * /oIPf 

1 — Pn 



(69) 



Pi = 



16ir 



, -2p T *l<r* r 2 / 1 - P(A . -p,»IU* J / 1 ~ Pl \ 

+ C T in I Pr ., ) + c h)[PF 1 — I 



70: 



2e (p r +p F )2 /2 ,2 /o . ^ 



/oP 



r2 /n 2 (p,^ Pn )/n(p, 



1 — P o 
rj 2 

1 — Pn 



Case 0: "P«7*i Off" 

We want to evaluate the integral 



Pa = J Pz( -I',//) t/.f <fy, 



(71) 
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where o 2 is the surface outside of a circle of radius pa with center at, the 
origin of coordinates, and pz(x,y) is obtained from (52) by making 

-4=0 (72) 

in the expression 

r = \/{x - AY + y 2 - 

After changing the variables according to (56), the \J/ dependence dis- 
appears from p 2 (x,y) and the probability of error (71) is 

P 2 = 2tt f th(r)r dr. (73) 

The most general term of this integral is proportional to 

k-j[^/.(5)/.(5)/.(5>*. m 

Over the range of integration, 

£»l; (75) 

also, for a small variation of r, the exponential varies much faster than 
the modified Bessel functions. Consequently, 



Vt S h I fiS ) /. ( fi? ) /. ( fi? ) / C -'""' r * 



-p »/2<r* j / PlPO \ y f P2P0 ^ j ( PiPl) 



(I 

Substituting this result in (73), we get 

P0 2/2,2 r- 



(70) 



1 + 2«-" ,/ * 1 h ( P ;f ') + if**'" 1 /o 2 ( P -^°) 
+ e -" 2/2 '' /„ (™j\ + 2 e - ( " + " )W /o fe) h fe°) (77; 
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